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We study the optimal diffusive transmission and absorption of broadband or polychromatic light 
in a disordered medium. By introducing matrices describing broadband transmission and reflection, 
we formulate an extremal eigenvalue problem where the optimal input wavefront is given by the cor¬ 
responding eigenvector. We show analytically that a single wavefront can exhibit strongly enhanced 
total transmission or total absorption across a bandwidth that is orders of magnitude broader than 
the spectral correlation width of the medium, due to long-range correlations in coherent diffusion. 

We find excellent agreement between the analytic theory and numerical simulations. 


One exciting development in optics in recent years is 
the coherent control of diffusing light in a disordered 
medium by shaping input wavefronts using a spatial light 
modulator (SLM) [1, 2]. Initially the emphasis was on 
using wavefront shaping (WFS) to focus light onto a 
wavelength-scale region (speckle) behind or within the 
disordered medium [3, 4], with potential applications for 
imaging; the optimal input wavefront in this case can be 
found by a simple sequential optimization of each pixel 
on the SLM, since each contributes to the local field at 
the focal spot independently. More recently there has 
been progress in the more challenging problem of opti¬ 
mizing global properties of the fields, such as the total 
transmitted power through the medium [5-7]. Motiva¬ 
tion came from theoretical concepts first formulated in 
the context of mesoscopic electron transport and localiza¬ 
tion theory [8-11], where it was predicted that in a loss¬ 
less diffusive medium there would always exist sample- 
specific “open channels” that will be transmitted almost 
perfectly. A closely related effect is the coherent enhance¬ 
ment of absorption (CEA) to near unity via WFS in a 
disordered medium that on average only absorbs a small 
fraction of the input light [12, 13]. Incomplete control 
of the input wavefronts reduces the possible enhance¬ 
ments [14, 15], but large enhancements are still observ¬ 
able under realistic conditions [6, 13]. 

The physical basis of these coherent control effects is 
manipulation of the multiple-scattering interference in 
the medium to violate the expected behavior for inco¬ 
herent diffusion. Hence these effects would seem to be 
intrinsically narrowband, limiting their applications in 
contexts such as power delivery, communications, or en¬ 
ergy conversion, in which larger bandwidths may be re¬ 
quired. The expected bandwidth is limited by the fre¬ 
quency correlation scale, Sco, which for lossless transmis¬ 
sion is the inverse of the time to diffuse across the thick¬ 
ness L of the medium, Soj « D j{D = Ic/d is the 
diffusion constant in d dimensions, and I is the transport 
mean free path) [16]; for CEA, duj « c/G, where la is the 
ballistic absorption length [12]. Eor a broadband signal 


with bandwidth Aw 3> Sco, a natural hypothesis is that 
the effective number of independent frequencies would 
be Mgff « 1 -|- Alo/Slo, and that the maximal achievable 
enhancement decreases as 1/Meff. Indeed this is exactly 
the behavior found in experiments maximizing the focal 
intensity of polychromatic light on a single speckle spot 
using SLMs [17-20]. However we will show that this is 
not the case for the total transmission or absorption due 
to the long-range spectral correlations of coherent diffu¬ 
sion [21-24], which are unimportant for speckle statis¬ 
tics but play a major role for the global properties [6]. 
These correlations dramatically reduce the effective num¬ 
ber of independent degrees of freedom, and instead of 
the linear scaling, we find Mgg « y/ Alo/6lu, allowing 
substantial coherent control of transmission and absorp¬ 
tion over large bandwidths. For example, for a lossless 
diffusive sample with average 2% transmission, the total 
transmission can be enhanced 10 times across bandwidth 
Aw « 60(lw; similarly for a thick diffusive sample with 
average 3% absorption, the total absorption can be en¬ 
hanced 10 times across Aw « 60(5w. 

We begin by defining a broadband flux matrix, based 
on the monochromatic transmission matrix t(w) that 
relates the incident field J/in) to the transmitted field 
l/t) = t|/in); the field vectors are written in the basis 
of N input and output modes carrying unit flux, and 
we assume N ^ 1. The monochromatic transmitted flux 
(V'tj/’t) is the expectation value (/inlt'i(w)t(w)lV;in) of the 
Hermitian matrix tH; it follows that the most open chan¬ 
nel for monochromatic light corresponds to the largest 
eigenvalue of tH [8-11, 25-30]. For polychromatic light, 
the role of tH is replaced by 

A = y (iw/(w)ti(w)t(w), (1) 

where /(w) is the power spectrum of the incident light 
normalized to / (iw/(w) = 1. When the transmitted 
flux is measured with a sufficiently-long integration time, 
beating between different frequencies averages away, and 
the total transmission for incident light with spectrum 
/(w) and wavefront \Hin) is simply (/injAj/in). Since A 
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is still Hermitian, the optimal wavefront is again given by 
the eigenvector with the largest eigenvalue. A broadband 
reflection flux matrix can be defined similarly, with r^r 
replacing tH. Note that a monochromatic open channel 
at some frequency within the spectral envelope is gener¬ 
ally no longer an eigenvector of A, so we immediately 
know that it will not provide the optimal broadband 
transmission. The optimization figure of merit here is 
essentially the frequency-averaged transmission with the 
bandwidth and weight specified by /(w); as we shall see 
later in Fig. 2(b), the optimal wavefront has its transmis¬ 
sion enhanced rather uniformly across the target band¬ 
width, so optimizations aiming for uniformity (such as 
maximin) will yield similar results. 

In the diffusive regime (A ^ A, where A is wave¬ 
length), each matrix has a bimodal eigenvalue 

density Pttt(T) = T/{2Ty/l — T) where T is the average 
transmission [8-11]. The distribution has support up to 
r = 1, meaning monochromatic open channels always 
exist for ^ 1 in the diffusive regime. Since the trans¬ 
mission matrices at different frequencies do not commute, 
the eigenvalue density of the broadband matrix A will be 
very different, as we show below. 

First, we study the simpler situation when A is given 
by a sum of matrices at discrete frequencies that we as¬ 
sume are so widely separated that correlations between 
them are negligible. Subsequently we will adapt this the¬ 
ory to treat a continuous input spectrum with long-range 
correlations, relevant for more experiments. Hence ini¬ 
tially we take I{lo) = X[m=i ^m<5(w — ojm) and assume 
no correlation between the M matrices, {t{uJm) = tm}- 
The setup for M = 2 is illustrated in Fig. 1(a). The 
eigenvalue density for a sum of large, mutually uncorre¬ 
lated, non-commuting random matrices can be treated by 
methods developed in free probability theory, which gen¬ 
eralizes the concept of statistical independence to such 
matrices [31, 32]. Specifically one can apply an addition 
rule [33] to find an implicit equation for the eigenvalue 
density of their sum. For the matrix A, define gA^z) as 
the Stieltjes transform (resolvent) of the eigenvalue den¬ 
sity PA'-, applying the addition rule, one finds that the 
unknown resolvent gA can be obtained from the follow¬ 
ing implicit equation (details in [34]) 

M — 1 ^ 

" + VJz) = ^ i^rngAiz)), ( 2 ) 

with the known resolvent gp ^ that is determined from 
the bimodal distribution ^ . We then apply standard 
root-finding algorithms to this equation to find gA{z) and 
obtain the desired eigenvalue density through the inverse 
Stieltjes transform pa{T) = — lim£_,,o+ Img^(T -|- ie)/7r. 
Results for the general (M > 2) cases are given in Fig. SI 
in [34]. Here we examine the simpler M = 2 case for dif¬ 
ferent combinations of weights {VFi, IF 2 }, shown as solid 
curves in Fig. 1(b) (here, Ti = 0.027, T 2 = 0.021). The 



FIG. 1: Total transmission through a disordered slab for inci¬ 
dent light with two discrete frequencies, (a) Schematic setup, 
with a disordered slab in a multimode waveguide and poly¬ 
chromatic light incident with a shared wavefront, (b) Density 
of the polychromatic transmission eigenvalues as calculated 
numerically by solving the wave equation (symbols) and an¬ 
alytically using Eq. (2) from free probability theory (lines). 
Numbers indicate the intensity weight Wi. 

Wi = 0, W 2 = 1 case corresponds to the monochro¬ 
matic bimodal distribution. With increasing Wi, the up¬ 
per edge Tniax decreases as expected; the residual peaks 
near T « Wi and T « W 2 can be traced back to the open 
channels of the constituent matrices Wit\ti and W 2 t\t 2 - 
From the case Wi = W 2 = 1/2, we see that Tmax ~ 0.59 
is larger than the (1 -I- T')/2 « 1/2 one would obtain 
from using the monochromatic open channels as the in¬ 
put wavefront. However, we note that in the limit of 
Ti,T 2 -a 0, Tmax will still approach 1/2. 

We perform numerical simulations to validate the an¬ 
alytic prediction. As illustrated in Fig. 1(a), we simulate 
a 2D disordered slab of thickness L and width IF = 3L in 
a waveguide geometry with background refractive index 
no = 1.5 and slab permittivity e(r) randomly sampled 
between Uq ± 0.9 at each grid point. Using the recur¬ 
sive Green’s function method [35], we obtain the N-hy-N 
transmission matrix (here N = 647) of the wave equation 
[V^ + (a;/c)^e(r)]'!/(r) = 0 for 600 realizations of disorder, 
at two frequencies wi = 390c/L and 0 J 2 = 410c/L (aver¬ 
age transmissions Tf = 0.027, T 2 = 0.021; the variation 
of N is negligible) that are much further than Suj apart 
(here UJ 2 — 1^1 ~ 290(5w). The resulting eigenvalue den¬ 
sities of the two-frequency matrix A, shown as symbols 
in Fig. 1(b), agree perfectly with the analytic prediction 
with no fitting parameters. 

To see the effect of transmission correlations, we per¬ 
form wave simulations of the diffusive medium for a 
broadband input with uniform spectral weights I{oj) 
over bandwidth Aw, centered at wq = 400c/L (where 
T = 0.025). The numerically obtained maximum eigen¬ 
value Tmax of the broadband matrix A is plotted as 
blue circles in Fig. 2(a), as a function of Aoj/Soj, where 
Suj = 0.069c/L « 21D/L^ is defined as the full width 
at half maximum (FWHM) of the transmission spec¬ 
trum for the monochromatic open channel (black line in 
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FIG. 2: Broadband transmission open channels, (a) Maxi¬ 
mal eigenvalue Tmax of the broadband flux matrix A and the 
enhancement t] = T^ax/r obtained from numerical simula¬ 
tions (blue circles) and analytic theory accounting for long- 
range correlation (blue line), showing the highest achievable 
frequency-integrated transmission across bandwidth Atj. The 
two lines below show the would-be maximal transmission 
(green, dashed) and the transmission of the monochromatic 
open channel (orange, dot-dashed) if there were 1 -I- Alj/Suj 
uncorrelated frequencies. Black dotted line indicates the aver¬ 
age transmission, (b) Transmission as a function of frequency 
when the input wavefront is fixed to the optimal eigenvector 
with different bandwidths Au. 


Fig. 2(b)); note that the FWHM of the open channel 
coincides with the FWHM of the speckle intensity corre¬ 
lation (see Fig. S2 in [34]). In Fig. 2(a), we find that at 
all bandwidths, T^ax (blue circles) is much larger than 
the prediction of the uncorrelated model when the effec¬ 
tive number of independent frequencies is taken as Meg = 
1 -I- Auj/Suj (green dashed line), which itself is larger than 
the frequency-averaged transmission of the monochro¬ 
matic open channel, (1 -|- TAuj/Suj)/{1 + Aoj/Suj), when 
one assumes MgS = 1 + Auj/Suj (orange dot-dashed 
line). The transmission spectra of the optimal broad¬ 
band eigenvectors cover the target bandwidth rather uni¬ 
formly, as shown in Fig. 2(b) for representative band- 
widths. Clearly spectral correlations beyond the FWHM 
are critical and should allow much greater coherent con¬ 
trol of broadband transmission. 

To account for these spectral correlations, we adopt an 
approach similar to the treatment of spatial correlations 
in Ref. [6]. We hypothesize that even in the presence 
of spectral correlations, the eigenvalue density can still 
be described by Eq. (2), but with M replaced by some 
effective number of independent frequencies, Mgg < 1 + 
Auj/Suj. We focus on the case where the spectral weights 
Wm is uniform, for which Eq. (2) takes a simpler form 


9a{z) 

MeS 


— 9tU 



MeS-l \ 
9a{z) ) ■ 


( 3 ) 


This coincides with Eq. (3) in Ref. [14] when A is taken 
to be iH with i being a “filtered” matrix that only has a 
fraction mi = l/Mgg of the input channels (columns) of 
the full matrix t, as in experiments where all output lights 
are measured but only a fraction of the incident channels 


FIG. 3: (a) Density of the broadband transmission eigen¬ 

values for various bandwidths, calculated numerically from 
simulations (symbols; bandwidths Auj/Suj — 0, 1.0, 2.4, 
5.9, 15, 40) and analytically from Eq. (3) with an effec¬ 
tive number Meft of independent frequencies (lines; Mefi = 
1.0,1.3,1.7, 2.4, 3.9, 6.6). (b) M^g as a function of the square- 
root bandwidth, evaluated numerically from Eq. (4) (symbols) 
and analytically from Eqs. (5)-(6) (line). 


is controlled by the SLM [6]. Given this equivalence, we 
can use a property of the filtered matrix [14] 

^ (4) 

Mes Var(r) ’ 

to determine Mgg, where r and t are the eigenvalues of 
A and of respectively. With the broadband eigen¬ 
value density from simulations (symbols in Fig. 3(a)), we 
confirm that Eq. (4) provides the correct value of MgS 
that, through Eq. (3), predicts analytical eigenvalue den¬ 
sities (lines in Fig. 3(a)) that agree well with the numer¬ 
ical data. Mgg obtained in this manner scales with the 
square root of the bandwidth (circles in Fig. 3(b)). 

The quantity Var(f) can be expressed in terms of the 
disorder average of certain products of four transmission 
amplitudes tab, and the disorder averages can be car¬ 
ried out analytically using impurity-averaged perturba¬ 
tion theory (details in [34] and Figs. S3-S4). We find 

1 _ f f dujiduj2 C^'^'>{uji,uj2) , . 

Mgg JJ Aw2 C(^)(a;o,a;o) ’ 

where {uji,uj 2 ) is the mean-normalized spectral cor¬ 
relation {Ta{uJi)Ta{uj 2 ))/{Ta{uji)){Ta{uj 2 )) -1 of the total 
transmission Ta — J2b with the brackets denot¬ 

ing average over disordered samples; the dependence on 
mode index a drops out due to the normalization. In our 
system, is well described by (see Fig. S4(a) in [34]) 


{UJI,UJ2) 


I 

m' 


2 sinh(x) 
X cosh(a;) 


sin(x) _ ^ 
cos(a;) 


( 6 ) 


where x = \/2\uJi — uj 2 \L'^/D. In Eq. (6), the first term 
is the long-range correlation that decays as \uji — 

(see Ref. [23]), while the second term is a finite-T cor¬ 
rection [34, 36]. Eqs. (5)-(6) provide an analytic expres¬ 
sion to calculate Mgg without free parameters and per¬ 
fectly agrees with the Mgg obtained from simulations. 
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as shown in Fig. 3(b). The blue solid line in Fig. 2(a) 
is calculated with this analytic expression of Mes, and 
it explains the much larger potential transmission en¬ 
hancement through WFS than expected from the uncor¬ 
related model. Specifically, when Aw falls in the regime 
1 <C \/ Aw/Jw l/T, the relevant values of are 

dominated by the |wi—tail in the long-range con¬ 
tribution, giving rise to the scaling of Meff ~ Aw/Jw 
and a parametrically larger Tmax- Note that Eqs. (5)-(6) 
show that Meff and T^ax depend only on the bandwidth 
Aw and the average transmission T. 

In addition to enhancing transmission, coherent super¬ 
positions of the input modes can also enhance absorption 
(CEA). Consider a thick diffusive scattering medium with 
A ^ Z <C y/Ua/d < L, where la is the ballistic absorp¬ 
tion length. As the thickness L is larger than the diffu¬ 
sive absorption length La = ^Jlla/d, the transmitted flux 
is exponentially small, so any light that is not reflected 
can be considered absorbed. As Z ^ la, most incident 
light is reflected before it propagates far enough to be 
absorbed, so the average absorption is low. However, 
there exist eigenchannels that can be nearly completely 
absorbed at one frequency when the number of input 
channels {i.e. degrees of freedom to be controlled) is large 
enough that N'^l/la ^ 1 [12, 37]. The minimum reflec¬ 
tion (corresponding to the maximum absorption) is the 
smallest eigenvalue of rV, and in the N oo limit the 
monochromatic eigenvalues follow a known bimodal dis¬ 
tribution, PririR) = 2a-\/(l — R)/{aR) — l/(7r(l — RY), 
where a = Z/Z^ 1 [38, 39]. For broadband light with 
spectrum /(w), we instead look for the eigenvalues of A 
as defined in Eq. (1) just with t(w) replaced by r(w). 

We perform numerical simulations for the geometry 
shown in the inset of Eig. 4(a) with thickness L and 
width W = 0.43L and with a weak uniform absorp¬ 
tion Im(e) = 3 X 10“® in the diffusive medium (corre¬ 
sponding to a = 2 X 10“'*, R = 0.97, and N — 323 
near wq = 1400c/L). Again, we consider broadband 
incident light with uniform spectral weights /(w) over 
bandwidth Aw, and numerically evaluate the reflection 
matrices r(w) and the eigenvalues of the broadband flux 
matrix A. The maximum absorption, 1 — Rmin, is plot¬ 
ted as blue circles in Eig. 4(a) as a function of Aw/^w, 
where (5w = 0.14c/L « 12c/Za is defined as the FWHM 
of the absorption spectrum for the monochromatic CEA 
channel (see Eig. S5(a) in [34]). Similar to the broadband 
lossless transmission, here we find the maximal absorp¬ 
tion to be much larger than the prediction if one were to 
ignore long-range spectral correlation (green dashed and 
orange dot-dashed lines). 

The density of broadband reflection eigenvalues is 
well described by Eq. (3) with tH replaced by Ar 
and with M^s given by Eq. (4) (see Fig. S5(b) 
in [34]), confirming the hypothesis that one can use 
an effective number of independent frequencies to de¬ 
scribe the broadband eigenvalue distribution. Analyt¬ 



FIG. 4: Broadband coherently enhanced absorption (CEA). 
(a) Maximal frequency-integrated absorption 1 —Amin and the 
enhancement rj = {1 — Amin)/(1 — R) obtained numerically 
(blue circles) and analytically (blue line) across bandwidth 
Aw. The two lines below show the would-be maximal absorp¬ 
tion (green, dashed) and the absorption of the monochromatic 
CEA mode (orange, dot-dashed) if there were 1 -b Auj/Soj un¬ 
correlated frequencies. Black dotted line indicates the average 
absorption. Inset shows schematic setup of the system; a re¬ 
flecting boundary on the right ensures that transmission is 
zero and absorption is 1 — A. (b) Mefr as a function of the 
square-root bandwidth, evaluated numerically from Eq. (4) 
(symbols) and analytically from Eqs. (5) and (7) (line). 


ically, Meff is again given by Eq. (5), just with 
replaced by the spectral correlation C^^\uji,uj 2 ) = 
(Aa (wi) Aa (W2 )) / (Ao (wi)) ( Aq (W2 )) -1 of the total reflec¬ 
tion Ra — J2b kfcoP) which in our system is well described 
by (see Eig. S4(b) in [34]) 


(wi, W2) 


1 - A 

AA(1 -b A) 


2 

1 + y 


(1-A) 


(7) 


where y = Reiyi -b fjwi — W 2 lZa/c. Here, the first term 
decays as [wi — W 2 l“*'^^ and is the long-range reflection 
correlation derived in Ref. [24], while the second term is a 
correction for finite 1 —A [34]. Eqs. (5) and (7) provide an 
analytic expression for MgS and is plotted as the red line 
in Fig. 4(b), with its prediction of the maximal absorp¬ 
tion plotted as the blue line in Fig 4(a). When the band¬ 
width Aw falls in the regime 1 ^ -^/Aw/Jw ^ 1/(1 — A), 
the is dominated by the jwi— W2l“*/^ tail in the long- 
range contribution, and Mes scales as a/ Au}/5u), giving 
rise to large potential enhancements of absorption similar 
to the lossless transmission case. 

Using a single spatial wavefront to control polychro¬ 
matic light does introduce a loss of control, on top of the 
incomplete channel control [14] present in previous WFS 
experiments with narrowband light. But long-range spec¬ 
tral correlations significantly reduce this loss, making 
the coherent control of short pulses and multi-frequency 
laser beams potentially feasible despite the rather narrow 
spectral correlation width, Suj, of the open channels and 
CEA channels. The formalism for uncorrelated matrices 
(Eq. (2)) can also treat spatially incoherent light with 
multiple uncorrelated transverse modes or unpolarized 
light with two independent polarizations. 
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Supplementary Material 


Eigenvalue Distribution for Sums of Uncorrelated 
Matrices 


X 2 are asymptotically free, their _R-transforms are addi¬ 
tive [1, 2, 5] 


In this section, we provide details for obtaining the eigen¬ 
value distribution for a sum of uncorrelated matrices, 
namely Eq. (2) in the main text and its extensions. For 
a Hermitian matrix X with eigenvalue distribution px, 
its Stieltjes transform (also called resolvent) gx is [1, 2] 

9 x{z) = r ^^dz' (S.l) 

J_oo 


for complex-valued z. For example, for a transmission 
matrix tm in the diffusive regime {X I L), the eigen¬ 
value distribution of follows [3] 


PtltJT) 


0 , 


sech^(2—) < T < 1 

otherwise 


(S.2) 


where Tm is the average transmission, so its Stieltjes 
transform is 




1 

2 


T 

— Arctanh 

zyl - 


Tanh(l/7;„) 

\/\- z 


(S.3) 


For the reflection matrix Vm in the diffusive regime with 
weak absorption that we consider in CEA (A ^ ^ 

\/lI^ < L ), the eigenvalue distribution of r^^rm follows [4] 


(R) 


/l_ R _ , n ^ 7 ? < ^ 

7ril-R)^y amR ’ - 1 + a^ 

0, otherwise 

(S.4) 


where = l/la = — Rm^/4:Rm, so its Stieltjes trans¬ 

form is 


- , ^ 1 + 2am atz , , 

9rlrJ^) = - 

Given the Stieltjes transform, we can apply Cauchy’s in¬ 
tegral formula to get back the eigenvalue distribution 

Px{z) = -- \im Imgxiz + ie). (S.6) 

TT e-> 0 + 

In general, knowing the eigenvalue distribution of two 
individual matrices is not enough to infer the eigen¬ 
value distribution of their sum, as the two matrices do 
not share the same eigenbasis. However, free probabil¬ 
ity theory identifies a sufficient condition, called asymp¬ 
totic freeness, under which the eigenvalue distribution of 
the summed matrix cab be computed. When Xi and 


Rxi+X2{y) = Rxiiy) + Rx^iv), (S.7) 

where the i?-transform is defined as 

Rx{y)=9x^{y)-^, (S.8) 

y 

with g^^ being the functional inverse of gx- 

For a discrete spectrum /(w) = ~ Wm), 

Eq. (1) of the main text becomes 


A = Y,Wmtltm, (S.9) 

m 


where tm = t(ujm)- In general, two Hermitian random 
matrices of size N are asymptotically free in the limit 
iV —)■ oo if the two matrices are uncorrelated to each 
other. This is the case for the set of matrices {Wmtlntm} 
in Eq. (S.9) if the frequencies {ujm} are much further 
than Suj apart. Applying Eq. (S.7) to Eq. (S.9) yields 


I 

9Aiz) 


E 


Wmgp\ {Wm9A{z)) 


1 

9a{z)_ ’ 


(S.IO) 

which is Eq. (2) in the main text. Given the weights 
{Wm} and the average transmissions {Tm} (or {Rm} for 
the reflection case), we can numerically solve Eq. (S.IO) 
for gA^z) using the known g^t ^ from Eq. (S.3) (or ^ 
from Eq. (S.5)), and obtain the eigenvalue distribution 
PA through Eq. (S.6). 

The minimal and the maximal eigenvalues also follow 
from Eq. (S.IO). At the lower or upper edge Zq of the 
distribution, the probability density becomes zero, and 
its slope diverges. That means gA^Zo) is purely real, and 
I ( 7 ( 4 ( 00 ) I = 00 . Taking the z derivative on both sides of 
Eq. (S.IO), we see that \g'A{zo)\ = 00 when 


1 _ _l_ I 

^'a(^o) “ ^ ^ iWmgAizo))) 9\{zo) 

(S.ll) 

We can solve Eq. (S.ll) for the real-valued gA{zo), and 
plug gA{zo) back into Eq. (S.IO) to obtain zq. In general, 
there will be two solutions of gA^zo), which gives the 
minimal and the maximal eigenvalues. 

The unknown in Eq. (S.ll) is real-valued and can be 
obtained with basic bracketing methods. When one edge 
Zq is found, we can start from z = Zq with the known 
gA{zo) being the initial guess, and use Newton’s method 
to solve Eq. (S.IO) for the complex-valued gA{z), with 
z changing incrementally based on the magnitude of the 
Jacobian. 
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FIG. SI: Distribution of transmission eigenvalues for a set of 
discrete and uncorrelated frequencies with weights following 
a discrete Gaussian distribution with standard deviation a. 
Symbols: simulation. Lines: solution of Eq. (S.IO). 



FIG. S2: Spectral correlation of transmission through disor¬ 
dered media, for the lossless system considered in Figs. 1-3 of 
the main text, (a) Transmission spectrum for the monochro¬ 
matic transmission eigenchannel with the highest transmis¬ 
sion at oJo- (b) Spectral correlation function of the speckle 
intensity. 


In Fig. SI, we plot the solution of Eq. (S.IO) with 
the weights {Wm} given by a discrete Gaussian distri¬ 
bution with standard deviation a. To validate the ana¬ 
lytical prediction, we carry out wave simulations of dif¬ 
fusive media (see main text for system parameters; here 
oj = 410c/L, T = 0.021), with the uncorrelated matrices 
given by the transmission matrices of different disorder 
realizations. The numerically calculated eigenvalue den¬ 
sities are shown as symbols in Fig. SI and agree perfectly 
with the analytic prediction with no fitting parameters. 

The effects of incomplete channel control, common in 
experiments, can be included in our formalism by replac¬ 
ing tm and Vm with the filtered matrices tm and fm in¬ 
troduced in Ref. [ 6 ]. 


Spectral Correlation Width for Transmission 

In Fig. S2(a), we show the transmission spectrum of 
the monochromatic open channel, for the lossless sys¬ 
tem considered in the main text; its full width at half 
maximum (FWHM) is what we define as the spectral 
correlation width 6uj. As shown in Fig. S2(b), this 
Soj coincides with the FWHM of the spectral corre¬ 
lation function of the speckle intensity, Cabi^iTi^ 2 ) = 
{Tabit^i)Tab{^^2)) / {Tahiti)) {Tab{^2)) “ 1, where Tab = 
ItfeaP- 


Relating Mgs to Spectral Correlation 

Here we derive and validate Eq. (5) in the main text 
and its more general form, which shows that the effective 
number Meg of uncorrelated frequencies is simply the 
inverse of the frequency-averaged spectral correlation of 
the total transmission or total reflection. To compute 


Meff = Var(T)/Var(f) as defined in Eq. (4) of the main 
text, we need to evaluate 


Var(f) = -(Tr(AA))- (-(Tr(A)) 


N 


jj du}iduj2l{oJl)I{ u}2)TiT 2F^^^ (uJi,UJ2) 


(S.12) 


where 


(a;i,a;2) 


(Tr(4ti4t2)) 

NT 1 T 2 


(S.I3) 


ti = Ti = (Tr(4ti))/iV, and the brackets denote 

averaging over different disordered samples. 

The quantity contains {Tr{t\tit\t 2 )) = 

Ea.h.a',b'(4a(wi)tt„,(a;i)4/a'(w2)tJ/a(w2))- Mean¬ 

while, the spectral-and-channel correlation of the 
speckle intensity involves a closely related quantity 
{Tab{.^l)Ta'b'{.<^2)) = (4a(^1 )tj^(wi )tf,'a' (^ 2 )^ 6 /^/(^ 2 )) , 

which in the diffusive regime separates into three 
contributions, Gi, G 2 , and G 3 [7, 8 ] (also see reviews 
in Refs. [9-12]). In waveguide geometries where the 
incident and outgoing channels are quantized into 
waveguide modes, the channel dependences appear as 
Kronecker delta symbols [ 8 , 13], and the Ci and C 2 con¬ 
tributions of {Tab{uJi)Ta>b> {UJ 2 )), shown diagrammatically 
in Fig. S3(a), give 


(4a(wi)tbQ(wi)4'a'(w2)t6/a'(‘*^2)) — {Tab) {Ta'b')'X 

[1 + 5aa'Shb'Ci{uJi,Ul2) d" (^aa' + 4b')42 (^1, W 2 )] ■ 

(S.14) 


The G 3 contribution is negligible for our purpose here. 
Also, we drop the weak frequency dependence of {Tab) 
and {Ta'b')- The four-amplitude average in the trace 
of Eq. (S.I3) differs from Eq. (S.14) only by switching 
the indices a and a' on the conjugated amplitudes; dia¬ 
grammatically this corresponds to bending the legs of the 
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FIG. S3: Diagrams showing the leading-order contributions for (a) the term a'{s^ 2 )tya,'['-^ 2 )) in the speckle 

intensity correlation, and (b) the term {tba(u>i)tbai {<-^i)tb'a' {^2)tb'ai‘^'2)} in ths trace in the definition of Red and blue colors 

indicate propagation at uji and UJ 2 respectively. Solid and dashed lines correspond to retarded Green’s function (unconjugated 
transmission coefficient) and advanced Green’s function (conjugated transmission coefficient), respectively. Vertical black lines 
indicate the boundary of the disordered medium. 


dashed lines on the incident side, as shown in Fig. S3(b), 
and the effect is to swap 5aa' and 1 in each term, giving 

= (Tab) {Ta'b') X 

[Saa' + Sbb'Cl (U! 1 , 012 ) + (1 + Saa'Sbb')C2{uJi,UJ2)] ■ 

(S.15) 

Summing over a,a',b,b', the fourth term becomes 
negligible, and we get F(^)(wi,a; 2 ) = Ci(wi,a; 2 ) + 
NC 2 {oji,oj 2 )- Instead of C 2 , it is easier to 
use the mean-normalized correlation { 001 , 0 ) 2 ) = 
{Ta{oJi)Ta{o) 2 ))/{Ta{ 0 Ji)){Ta{ 0 J 2 )) - 1 of the total trans¬ 
mission Ta = Setting a = a' and summing 

over b and b' in Eq. (S.14), we see that C'C") { 001 , 002 ) = 
{l/N)Ci{ooi,oo 2 ) + C' 2 (a;i,a; 2 ), so we get 

F('^)(u;i,a;2) = (S.16) 


To test the validity of Eq. (S.16), in Eig. S4(a) we plot 
fC") (blue circles) and NC^C (black crosses) calculated 
from numerically obtained transmission matrices of the 
lossless system considered in the main text. The agree¬ 
ment remains excellent even when absorption is included 
(data not shown). 

Inserting Eq. (S.16) into Eq. (S.12) and dropping the 
weak frequency dependence of T, we get 


Var(f) 

Var(T) 


JJ doJidoj2l{ooi)I{oj2) 


{001,002) 
CC") (wo, Wo) 


(S.17) 


Eq. (S.17) holds for arbitrary spectrum I{ 00 ) as long as 
T is roughly constant within the bandwidth and that 
the system is in the diffusive regime. However, the re¬ 
lation Mgff = Var(r)/Var(r) is only meaningful when 
the eigenvalue density of the broadband matrix A fol¬ 
lows the equal-weight uncorrelated model (Eq. (3) in the 
main text), which is not necessarily the case when the 
spectrum I{ 00 ) exhibits complex structures. When I{ 00 ) 




FIG. S4: Numerical validation of and of 

the analytical expressions for for transmission in the 

lossless system considered in Figs. 1-3 of the main text (a) 
and reflection in the weakly absorbing system considered in 
Fig. 4 of the main text (b). Blue circles are numerically cal¬ 
culated F^C and F^^^ as defined in Eq. (S.13). Black crosses 
are numerically calculated spectral correlation of the total 
transmission and total reflection, multiplied by the number 
of modes N. Red solid lines are the analytical expressions in 
Eq. (S.21) and Eq. (S.23). 


is a smooth function without large variation, the equal- 
weight uncorrelated model works well, and we can use 
Eq. (S.17) to obtain Meg. Eq. (5) in the main text is the 
special case when I{ 00 ) is a uniform distribution across 
bandwidth Aw. 

Similarly, for the reflection matrices, we have 

F(^)(wi,W2) = NC^^\oji,oo2), (S.18) 

where F(^^(wi,W 2 ) = (Tr(r|rir|r 2 ))/(A.Ri.R 2 ) — 1) and 
(Wi, W2) = (i?Q(wi)i?a(w 2 ))/(i?o(wi))(i?Q(w 2 )) — 1 
with the total reflection being Ra — J2b kbaP- Although 
the reflection correlation is more complicated involving 
the effect of coherent back scattering [14, 15] and of short 
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optical paths [16], these effects are present in both 
and so the final expression Eq. (S.18) appears to 

be unchanged. Indeed, in Fig. S4(b) we numerically val¬ 
idate Eq. (S.18) for the weakly absorbing system consid¬ 
ered in the main text, and perfect agreement is observed. 
Eq. (S.18) also holds in the absence of absorption (data 
not shown). Therefore, Eq. (5) in the main text is also 
valid for reflection, just with replaced by 

In the absence of absorption, we can derive Eqs. (S.18) 
from (S.16) or the other way around based only on energy 
conservation. Energy conservation requires that rV -I- 
tH = I at both frequencies, so {TT{t\titlt 2 ))/N — T 1 T 2 = 
(TT{r\rirlr 2 ))/N-R 1 R 2 . Meanwhile, .Ro(w)-l-To(w) = 1 
for every incident channel a, so TiT 2 C^'^\(jJi,(jJ 2 ) = 
, a; 2 ). These two requirements demand 
Eq. (S.18) to follow from Eq. (S.16). 

When LOi = UJ 2 , the trace in is simply the sec¬ 

ond moment of the eigenvalues, which we can evaluate 
directly from Eq. (S.2) and Eq. (S.4). In such case, 
Eqs. (S.16)-(S.18) yields 


(wo, Wo) 


1 

Wf 



C^^'> (wo,wo) 


1 - R 
N{l + R)’ 


(S.19) 


Analytical Expression for Transmission 


Refs. [19, 20] derived the spectral correlation of the total 
transmission, as 


(wi, W2) 


1 2 sinh(a;) — sin(a;) 
NT X cosh(x) — cos(a;) ’ 


(S.20) 


where x = ^2|a;i — U! 2 \L‘^/D, L is the sample thickness, 
and D is the diffusion coefficient. For large frequency dif¬ 
ference, has a long-range |wi — 072 !”^^^ dependence. 
Various corrections has been derived to account for the 
effects of absorption [21, 22], finite illumination area [23], 
and interfacial layer [24, 25]. However, Eq. (S.20) pre¬ 
dicts that = (2/3)/(VT), which differs from 

Eq. (S.19) by a shift that vanishes with T. Indeed, 
from the numerically calculated transmission matrices, 
we find that Eq. (S.20) differs from the numerical results 
by a constant shift that vanishes when T —>■ 0. Such a 
constant shift was observed experimentally in Fig. 5 of 
Ref. [23], although it was treated as an artifact of the av¬ 
eraging procedure. Empirically, we find that the follow¬ 
ing expression accounts for the constant shift and yields 
the correct value of C^'^\ujo,ujo), 


{UJI,UJ2) 


1 

Wr 


2 sinh(x) 


X cosh(x) 


sin(a:) _ ^ 
cos(x) 


(S.21) 


for the systems of interest here. We will use these equal¬ 
ities for the analytical expressions of and in 
the next two sections. 

We also note that when oji = 0 J 2 , Eqs. (S.16) reduces 
to a simple statement Var(T) = War(To) that is a di¬ 
rect consequence of the randomness of the matrices and 
the small fluctuations of the eigenvalues with respect to 
disorder realizations. The transmission matrix t has sin¬ 
gular value decomposition t = where r is a 

diagonal matrix of the transmission eigenvalues, and U, 
V are unitary matrices. Then, (T^) = 
i^a) = The eigenvalues {r„} 

are not correlated with the matrix V, and the sum of 
eigenvalues (NT) has a fluctuation with respect to dis¬ 
order that is of order unity [17] and therefore negligi¬ 
ble when NT is large, so it suffices to take the ensem¬ 
ble average over the unitary matrix V. We further in¬ 
voke the isotropic approximation that all incident chan¬ 
nels a are equivalent, so the ensemble average becomes 
the average over the invariant measure of the unitary 
group, for which (|VanP) = ^/N and (|VonP|KmP) = 
(1 + Snm)/iN{N + 1)) {e.g., see Refs. [2, 18]). This 
gives Var(T,) = (T^) - = ((r^) - f^)/{N + 1) = 

Var(T)/(V -I- 1). Note that here Var(To) is the variation 
with respect to disorder realizations (while a is fixed), 
whereas Var(T) is the variation with respect to the eigen¬ 
value distribution. This derivation also applies to the 
reflection matrices. 


This is Eq. (6) in the main text. This expression is plot¬ 
ted as the red solid line in Fig. S4(a). The only unknown 
parameter in Eq. (S.21) is the diffusion coefficient, which 
we set to D = 3.3 x 10“^cL. Once D is known from the 
the transmission correlation, we can use Eq. (5) in the 
main text to obtain analytically. 

Analytical Expression for Reflection 

Rogozkin and Cherkasov [16] derived the spectral corre¬ 
lation of the total reflection for a scattering medium with 
weak absorption (a = \/l/la —> 0). Taking the L » \/U^ 
limit of Eqs. (B5a) and (B6) in Ref. [16], we obtain 

where y = Re-\/l -I- i\uJi — W2|^a/c. As in the transmis¬ 
sion case, Eq. (S.22) exhibits long-range \u>i — 
dependence for large frequency difference. However, 
Eq. (S.22) was derived for very weak absorption (a —> 0, 
and R —)■ 1), and it only agrees with Eq. (S.19) in this 
limit. Indeed, from the numerically calculated reflection 
matrices, we find Eq. (S.22) to differ from the numerical 
results by a constant shift that vanishes when R —)■ 1, 
similar to the transmission case. Empirically, we find 
that the following expression accounts for the constant 
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FIG. S5: (a) Absorption spectrum for the monochromatic 

reflection eigenchannel with the highest absorption at ljq. (b) 
Density of polychromatic reflection eigenvalues for various 
bandwidths, calculated numerically from simulations (sym¬ 
bols; bandwidths Auj/Su — 0, 1.2, 2.7, 6.6, 17, 45) and 
analytically from Eq. (3) of the main text with an effec¬ 
tive number Mes of uncorrelated frequencies (lines; MeS = 
1.0,1.3,1.6,2.0,2.9,4.4). 


shift and yields the correct value of C'^'^)(a;o, wo)) 


^2^^ (Wl, W 2 ) 


1 - R 

N{1 + R)R 


2 

l + y 


(l-K) 


(S.23) 


This is Eq. (7) in the main text. This expression is plot¬ 
ted as the red solid line in Fig. S4(b). The only unknown 
parameter in Eq. (S.23) is the ballistic absorption length, 
which we set to la = 83L. Once la is known from the the 
reflection correlation, we can use Eq. (5) in the main text 
to obtain Mes analytically. 


Broadband Coherently Enhanced Absorption 

Here we provide additional data for the weakly absorbing 
system considered in Fig. 4 of the main text for broad¬ 
band CEA. Fig. S5(a) shows the absorption spectrum for 
the monochromatic reflection eigenchannel with the high¬ 
est absorption (lowest reflection); its FWHM is what we 
use to define 6uj. In Fig. S5(b), we plot the distribution of 
broadband reflection eigenvalues for various bandwidths 
Aw, comparing the simulation results (symbols) and the 
prediction from the uncorrelated model (lines) with an 
effective number of uncorrelated frequencies. 
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